id dl 




13a iSI 


fl/W- /0 5'78 



SIMPLE METHOD FOR THE DETERMINATION OF EFFECTIVE RAIN CELL 
DIMENSIONS AND ORIENTATION 

D. B. Hodge 


Tin© Ohio State Umv®ir§Sfiy 



Bepordmentf ®{ ISlosJs’ksi Enfjinesring 
Columlbws, Obi® 43212 


TECHNICAL REPORT 2374-15 
September 1973 

Grant Number NGR 36-008-080 


CASE FILE 
COPY 


National Aeronautics and Space Administration 
Office of Grants and Research Contracts 
Washington, D.C. 20546 




NOTICES 


When Government drawings, specifications, or other data are 
used for any purpose other than in connection with a definitely 
related Government procurement operation, the United States 
Government thereby incurs no responsibility nor any obligation 
whatsoever, and the fact that the Government may have formulated, 
furnished, or in any way supplied the said drawings, specifications, 
or other data, is not to be regarded by implication or otherwise as 
in any manner licensing the holder or any other person or corporation, 
or conveying any rights or permission to manufacture, use, or sell 
any patented invention that may in any way be related thereto. 



A SIMPLE METHOD FOR THE DETERMINATION OF EFFECTIVE RAIN CELL 
DIMENSIONS AND ORIENTATION 

D. B. Hodge 


TECHNICAL REPORT 2374-15 
September 1973 

Grant Number NGR 36-008-080 


■ National Aeronautics and Space Administration 
Office of Grants and Research Contracts 
Washington, D.C. 20546 



ABSTRACT 


A simple method is proposed for the determination of effective rain 
cell dimensions and orientation. Two rain cell models have been con- 
sidered: the circular cell and the elliptical cell. In both models it 

is assured that the rain rate is constant throughout the cell, that all 
cell locations are equally likely, and that the cell dimensions depend 
upon the rain rate. It is then shown that the effective cell dimensions 
and orientation may be deduced from rain rate statistics accumulated at 
two or three closely spaced sites, depending upon the model chosen. The 
results of this study along with previous estimates of rain cell dimen- 
sions indicate that the rain rate measurement points should be spaced on 
the order of one-half kilometer apart. The proposed method is simple 
and relatively inexpensive; thus, this approach is readily suited to the 
study of the dependence of rain cell size and orientation upon climatic 
region. The resulting rain cell characteristics are of direct value in 
the prediction of millimeter wavelength attenuation statistics on both 
single-terminal and diversity earth-space propagation paths as well as 
point-to-point terrestrial links. 
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I. INTRODUCTION 

The development of wi de-bandwidth communication systems has led to 
an increased demand for utilization of the frequency spectrum above 
10 GHz. Unfortunately, signals in this portion of the spectrum are 
attenuated by liquid water; and, in fact, severe attenuation may be en- 
countered when the signal passes through regions of intense rain rates 
such as those associated with severe thunderstorms ! 1 |. Thus, the 
overall communication system reliability may be ultimately determined by 
the spatial extent of those regions of intense rain rate through which 
the signal must pass and the probability of their occurrence. 

A simple, cylindrical rain cell model was studied in an effort to 
determine its utility in the prediction of attenuation statistics. Al- 
though this model grossly oversimplifies the physical characteristics of 
the actual rain cells, it was found to be quite useful for the prediction 
of long term attenuation statistics within one or two decibels' [2 j. This 
model has been generalized to an ellipsoidal shape in order to more 
accurately reflect the effect of cell shape on space diversity communi- 
cation systems where two spatially separated, parallel propagation paths 
are used to improve the overall system reliability. In the course of 
both of these modeling efforts it has been found that little data re- 
lating rain cell size, shape, orientation, and rain rate exists. Further, 
the dependence of such data upon climatic region is virtually nonexistent. 

The use of radar for the determination of these characteristics .'has 
been proposed and is being pursued. However, this approach is expensive, 
and it is unlikely that the instrumentation will be widely implemented 
in the near future. Therefore, an alternative approach for the measure- 
ment of these cell characteristics is proposed herein. This approach 
would utilize either two or three rain rate gauges spaced relatively 
closely, on the order of one-half kilometer, and a means of simultaneously 
recording the measured rain rates. Thus, the instrumentation require- 
ments would be minimal and the approach would be quite inexpensive to 
implement. 

The radar approach yields size, shape, and orientation information 
concerning each rain cell observed. This information must then be pro- 
cessed to obtain statistical cell properties. In contrast, the rain 
gauge approach will not yield characteristics of individual cells but 
will considerably simplify the data processing required to obtain the 
desired effective cell characteristics. The radar approach can also 
yield information concerning the vertical distribution of rain rates; 
whereas, the rain gauge approach yields only surface rain rate statistics. 
Initially, both techniques should be compared at one site in order to 
determine the degree of agreement or disagreement. 




II. THE RAIN CELL MODELS 

The rain cell models to be considered were chosen because of their 
tractabi 1 i ty . The intent was to represent rain cells in a phenomenal 
sense rather than to attempt to model individual rain cells. The first 
model considered was simply circular cylinder having a finite height 
above the earth's surface; the second model was an ellipsoid having two 
axes on the surface of a flat earth. The first model has two geometrical 
parameters, the diameter and the height, while the second model has four, 
the length, width, height, and orientation angle. Since the rain gauge 
method discussed herein deals only with surface rain rates, the vertical 
distributions of rain rate will be ignored. Hence, in the following, 
only the circular or elliptical surface cross sections of the cells will 
be considered. These geometries are shown in Fig. 1 where the x- and y- 
coordinate axes lie in the plane of the earth's surface. The orienta- 
tion of this coordinate system is rather arbitrary; but, for the sake of 




CELL GEOMETRIES 

Fig. lV The geometries of the circular and the 
elliptical rain cell models. 
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completeness, let us define the positive y-axis such that it is oriented 
in a northerly direction and the positive z-axis such that it is ori- 
ented upward. In both cases the coordinates X,Y denote the position of 
the center of the cell, and the remaining parameters are: 

d .= the diameter of the circular cell model, 

a,b = the major and minor axes of the elliptical cell model, and 

4 > = the orientation angle of the elliptical cell model with 

respect to the x-axis. *. 

The following basic assumptions are made in both cases: 

(1) The rain rate is constant throughout each cell. 

(2) In a statistical sense, there is a unique, single 
. valued relationship between the cell rain rate 

and the cell geometrical parameters. 

(3) All cell locations are equally likely. 

(4) Only one cell influences the rain rate measurements 
at any instant of time. 

Again, it must be stressed that the rain cells are to be modeled in a 
phenomenal, rather tha'n physical, sense. Note, also, that assumption (2) 
implies that d, a, b, and <j> are all dependent upon the rain rate, R. 
Finally, assumption (4) is in general too restrictive; however, in the 
cases to be examined here, the gauge spacings will be small compared to 
typical cell sizes. Therefore, since rain cells do not tend to occur 
immediately adjacent to one another, it is not likely that two different 
cells will simultaneously produce _rain at two closely spaced gauges. 

III. CIRCULAR CELL MODEL 

If a rain gauge is placed at the origin of the coordinate system, a 
record of the rain rate, R, may be obtained as a function of time. The 
percentage of time during which the rain rate is less than or equal to a 
given value, R Q , may then be readily determined from such data as a func- 
tion of R 0 - It, then, the data gathering period is sufficiently long, 
it will be assumed that this percentage distribution is equivalent to. 
the rain rate probability distribution. This distribution will be re r 
ferred to as the point rain rate distribution, P p +(R), where the rain 
rate subscript has been dropped for simplicity. p bince it has been 
assumed that all rain cell locations are equally likely, it follows that 
the point rain rate distribution does not depend upon the coordinate 
location of the measurement point. The point rain rate probability 
density function may then be defined in the usual manner, 

(1) Ppt( R ) = 

v d R 

The preceding pragmatic line of reasoning applies to all of the proba- 
bility density functions discussed in the following and, hence, will not 
be repeated. 

3 




A cell probability density function, p C ell( R )> will be associated 
with cells having rain rate R which are centered at X,Y. Again, Pcell(R) 
does not depend upon X or Y since the cells were assumed to be uniformly 
distributed. It should be reiterated at this point that a unique effec- 
tive cell diameter is assumed to be associated with each rain rate. 
Consequently, a cell of rain rate R centered within a circle of diameter 
d about the gauge wi 11 produce a rain rate R in the gauge, as shown in 
Fig. 2. This area about the gauge within which the occurrence of rain 



Fig. 2. Area of influence. 


cell center will produce rain in the gauge will be referred to as the area 
of influence, A pt (d). This area of influence depends upon the cell dia- 
meter and, thus, implicitly upon the rain rate. In fact, it can be noted 
that the size and shape of the cell and the area of influence are ident- 
ical for any given rain rate. Now, since the cells are uniformly 
distributed, it follows that the point rain rate and cell probability 
density functions are related by 


(2) Ppt( R ) “ A pt( d ) Pcell( R ) 
or 

(3) Ppt(R) - PcellW • 

If, now, a second gauge is placed a distance £ along the x-axis from 
the first gauge, where £ is less than the cell diameter, d, the geometry 
shown in Fig. 3 results. For any given rain rate the areas of influence. 
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Fig. 3. Two gauge geometry. 


A pt and A . , are identical. The joint area of influence. A- 
that area in which the occurrence of a rain cell center will p 
rain at both gauges simultaneously. From geometrical considerations, 
the joint area of influence is 

( 4 ) A joint( d »*) = sin" 1 jl- ^ . 

The joint rain rate probability density function, p- Q1 - t (R,ji) associated 
with the two gauges may now be related to the cell probability density 
function by 


p joint( R ’^ = A joint( d>£ ) Pcell^ 


oint 5 is 
roduce 


This joint rain rate probability density function is dependent upon only 
one rain rate since the model cons i de red permi ts the existence of only 
one rain cell at any instant of time. Thus, the simultaneous occurrence 
of dissimilar rain rates at the two gauges is not possible within the 
confines of this model. 

Substituting (4) into (5) yields 


(6) p joint (R ’^ = h ^ 2 sin " 1 / 1_ - £/d2 * 2 ' IPcellM 
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Letting 

< 7 > '*■ j 

and eliminating P ce n ( R) from Eq. (3) and (5) gives 

( 8 ) . . 2 WI=F] . 

Ppt< R ) 

Here the ratio of the two rain rate probability density functions is a 
relatively simple function of the geometrical parameters % and d. This 
functional relationship is shown in Fig. 4. Thus, if the two rain rate 
probability density functions are evaluated experimentally for a given 
value of rain rate and, of course, the gauge spacing is known, the 
effective rain cell diameter may be taken directly from Fig. 4. There- 
fore, the effective rain cell diameter, d(R), may be determined directly 
and simply from rain rate measurements at two sites. 

It is evident from Fig. 4 that the useful range measurement of y 
extends from roughly 1 to 20, i.e., d ranges from £ to 20£. The impli- 
cation here is that the gauge spacing should be small compared to the 
cell diameters of interest. Previous estimates of cell diameters indi- 
cate that they may range from 1 to 10 km with the smaller diameters being 
associated with the more intense rain rates [3] . Consequently, one may 
conclude that a gauge spacing of 0.5 km would be desirable for the im- 
plementation of this approach. Initially, this small spacing appears 
to be contradictory with rain gauge spacings cormionly used in the past. 
However, these measurements have generally been concerned with 
distributions of point rain rates rather than the interpretation of the 
data in terms of areal rain rate distributions. This distinction between 
point and areal rain rate statistics has been pointed out by Epstein ] 4 | . 

IV. ELLIPTICAL CELL MODEL 

The treatment of the elliptical cell model follows that of the * 
circular cell model with the exception of the increased number of degrees 
of freedom. The area of influence associated with a single gauge located 
at the origin of the coordinate system is shown in Fig. 5 and, again, is 
identical in size, shape, and orientation to the rain cell for any given 
rain rate. Locating an additional gauge at a separation distance of i 
along both the positive x- and y-axes yields the geometry shown in 
Fig. 6. The areas of influence associated with the gauges are again 
equal for any given rain rate, 

(9) A pt (a,b,<j>) = A r (a,b,<j>) = A 2 (a,b,<j>) = A 3 (a,b, <}>) . 
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RATIO OF RAIN RATE PROBABILITY DISTRIBUTION 
FUNCTIONS vs. GEOMETRICAL PARAMETER 
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Fig. 4. Ratio of rain rate probability distribution 
functions versus geometrical parameter y. 
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Fig. 5. Elliptical rain cell area of influence. 

And, the point rain rate probability density function is given by 
(10) Ppt( R ) = A pt (a,b,(|,) P ce11 (R) 

(U) P pt (R) = Ja b P cell (R) . 

The joint areas of influence are denoted by A.- (a,b,<j>,£') and are given 
by J 


ab 


(12) A 12 (a, b, 4 >,z) = — sin 


,-i 
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fcos 2 <f> sin 2 <|>\| 
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- 2 ..'(b 2 cos 2 6 + a 2 sin 2 <j>) 
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3 GAUGE GEOMETRY 

Fig. 6. 3 gauge geometry. 



(sin<|>+coS()>) 2 \l 


The ij subscript notation here and in the following obviously implies 
only the pairs: 1,2; 2,3; or 3,1. The three joint rain rate probability 

density functions may now be expressed as 


(15) Pj = A-j j(a,b,(j) ,«,) P ce n (R) - 
Eliminating p ce n(R) from Eqs. (11) and (15) yields 



(16) 


Pi j (R>£) 
Ppt^ R ) 




where the geometrical parameters are given by 



Comparing Eqs. (8) and (16), it is seen that the ratios of the proba- 
( bility density functions have the same functional dependence on the 

geometrical parameters for both models. Thus, given a set of joint and 
point rain rate probability density functions, the geometrical parameters, 
T j j , may be determined from Fig. 4. Eqs. (17), (18), and (19) may now be 
solved simultaneously to find a, b, and^in terms of this set of geo- 
metrical parameters, Y-jj: 






Note that the order of solution here is important since Eqs. (17) - ( 19 ) 
were not linear. The orientation angle, <j>, must be determined first from 
Eq. (20) and then used in Eqs. (21) and (22) to find the ellipse axes. 

The orientation angle, <)>, ranges from— it/ 2 to n/2; nevertheless, Eq. (20) 
still admits two possible solutions. This does not lead to an ambiguity 
since the value of a determined from these equations gives the length of 
the axis along the particular orientation angle chosen. _ Thus, a is hot 
necessarily the major axis of the ellipse, it is the axis oriented in 
the direction of the orientation angle, <j> , chosen from Eq. (20). 

Thus, given a set of joint and point rain rate probability density 
functions and the gauge spacing, it is possible to deduce unique 
effective ellipse axes and the ellipse orientation angle for each rain 
rate. The conclusions regarding practical gauge spacings discussed in 
the context of the circular model apply in this case as well. 

V. SUMMARY i 

A simple rain gauge method for the determination of effective rain 
cell dimensions and orientation has been described. Both circular and^ 
el 1 i pti cal rai n cell models have been treated. Either two or three rain 
gauges are required depending upon the model chosen. Further, the gauges 
should be spaced relatively closely, on the order of one-half kilometer. 

And, finally, the data processing required to extract the desired infor- 
mation from the rain rate statistics is quite simple and straightforward. 
Thus, a relatively simple, inexpensive technique is available for the 
statistical determination of effective rain cell dimensions and orientation. 
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